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The phase diagram of the massive chiral Gross-Neveu model (the 1+1 -dimensional Nambu-Jona- 
Lasinio model at large A'^) is investigated in the vicinity of the tricritical point. Using the derivative 
expansion, the grand canonical potential is cast into the form of a Ginzburg-Landau effective ac- 
tion. Minimization of this action by variational and numerical methods reveals both 1st and 2nd 
order phase transitions to a chiral crystal phase, separated by a tricritical line. These findings are 
contrasted to the massive Gross-Neveu model with discrete chiral symmetry where only 2nd order 
transitions have been observed. 
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I. INTRODUCTION 



Following Gross and Neveu 
fermion interaction models in 1 
grangians 



[H, we consider four- 
f 1 dimensions with La- 



C = iiJj^'ip — moTp'ijj- 



<r_ 

2 



(^V)2 



(GN) 
(NJL) 



(1) 



where we suppress flavor indices (i = 1...N). Only the 
large N limit is of interest to us, so that the classic 
no-go theorems 0, Q are bypassed. For vanishing bare 
mass mo — 0, the first variant (referred to as Gross- 
Neveu (GN) model hereafter) has a discrete chiral sym- 
metry ij: — > 75'(/', the second one (the two-dimensional 
Nambu-Jona-Lasinio (NJL) model Q) a continuous chi- 
ral symmetry ij: — > exp(i75Q!)^. The respective symme- 
try groups Z2 and U(l) manifest themselves in the bound 
state spectrum. Thus the GN model possesses massive 
kink and kink-antikink baryons, where the kink reflects 
most directly the underlying Z2 symmetry [5|. The NJL 
model features a pseudoscalar meson and a baryon both 
of which become massless in the chiral limit due to chiral 
symmetry breaking [1,0|- Whereas the meson ("pion") is 
a small fluctuation in the flat direction of the vacuum fa- 
miliar from the Goldstone theorem in higher dimensions, 
the massless baryon is a more exotic object specific for 
dimensions. It may be associated with a full turn 
around the chiral circle. Since this happens "adiabati- 
cally" as a function of a:, it does not cost any extra en- 
ergy and the baryon mass vanishes. One can show that 
the winding number of the chiral field is baryon number 
@. Hence the NJL baryon has a topological structure, 
reminiscent of the Skyrme model in 3-1-1 dimensions 

In the present work, we study the 1-1-1 dimensional 
NJL model at finite temperature T and chemical poten- 
tial /i. The main issue is the behavior of the chiral con- 
densates 



S 



mo 



P = V(^i75V'), (2) 



now regarded as thermal expectation values, as a func- 
tion of T and /x. If one enforces translational invariance 
and considers only homogeneous condensates as was done 



originally, the GN and NJL models appear to have iden- 
tical phase diagrams d, [l3| • From the physics point of 
view, this is puzzling, given that the difference in chiral 
symmetries has such a dramatic impact on the structure 
and mass of individual baryons. It was realized during 
the last few years that GN models at finite chemical po- 
tential may develop an x-dependent condensate, leading 
to solitonic crystal phases (for a recent review, see [Tlj). 
If one takes into account these crystal phases, the phase 
diagrams for the models with discrete and continuous 
chiral symmetry become quite different and the puzzle 
disappears. 

So far, the phase diagrams of the massless and mas- 
sive GN models have been worked out in full detail, to a 
large extent analytically [H, [l^. By contrast, thermo- 
dynamics of the NJL model including the crystal phase 
has only been explored in the chiral limit [1^ where it 
leads to the phase diagram shown in Fig. 1. Here, the 
horizontal line T ^ = e^/n (C « 0.5772 is the Euler 
constant) represents a 2nd order critical line in the (/i, T) 
diagram. Above this line, chiral symmetry is restored 
and the system behaves like an ideal massless Fermi gas. 
Below the line, chiral symmetry and translational invari- 
ance are broken. The condensate exhibits a helical shape 
(dubbed "chiral spiral" in Q), 



5- iP = me^''"^. 



(3) 



with a radius m depending only on T and a pitch e qual 
to 7r//i. This phase diagram has been derived in [l^ 
by "gauging away" the chemical potential through a lo- 
cal 75-transformation. Since this transforms the helical 
condensate 1^ back onto the standard mass m, the fi- 
independent phase boundary shown in Fig. 1 follows at 
once. Unfortunately, this trick does not work at finite 
bare fermion mass mo, and nothing is known yet about 
the phase structure of the massive NJL2 model. From 
Fig. 1 it is not obvious how the phase diagram will evolve 
as we turn on the bare mass. 

As a first step towards clarifying the phase structure 
of the massive NJL model, wc propose to zoom in onto 
the tricritical point and its neighborhood. The tricritical 
point is the point Ft at (^ = 0,T = Tc) where the hor- 
izontal critical line crosses the vertical line fi = 0. The 



2 



m = 



we minimize this action approximately by means of a 
simple variational calculation to gain some analytical in- 
sight. In Sect. IV we construct the phase diagram in the 
vicinity of Pt by means of a full numerical solution of the 
Euler-Lagrange equation. Sect. V contains a compari- 
son between the tricritical behavior of the GN and NJL 
models and our conclusions. 



chiral crystal 



II. GINZBURG-LANDAU THEORY NEAR THE 
TRICRITICAL POINT VIA THE DERIVATIVE 
EXPANSION 



1.5 



FIG. 1: Phase diagram of the massless NJL2 model in the 
{fi, T) plane according to Ref. [l4|. In the chiral crystal phase, 
the condensate has helical symmetry, see Eq. ([5]). 



segment of this latter line between T — and T = 
should also be viewed as a critical line. It lies at the 
boundary of the phase diagram just because the baryons 
are massless (the critical chemical potential at T = co- 
incides with the baryon mass). The identification of Pt 
as tricritical point will be fully justified in the course of 
our calculations where we shall construct a tricritical line 
ending at this point. 

Finite mo introduces a new axis into the phase diagram 
which we label by 



7 



(4) 



following earlier work on the GN model [Tl|, [l^ . This is 
a dimensionless, renormalization group invariant param- 
eter which replaces the bare fermion mass mo in physical 
observables. If we can understand the vicinity of the 
point Pt in (7, /i, T) space, we have a good chance to get 
basic information about the phase diagram such as or- 
ders of phase transitions and characteristics of different 
phases. On the other hand, the calculations in this region 
are significantly easier because the condensates are weak 
and slowly varying. This will enable us to avoid the te- 
dious numerical Hartree-Fock (HF) procedure and solve 
instead a Ginzburg-Landau (GL) type theory where the 
fermions have been "integrated out" , leaving behind an 
effective action for one complex boson field. Although 
such a calculation is approximate, it will become bet- 
ter and better as we approach the tricritical point. In 
particular it will enable us to determine rigorously the 
asymptotics of how the phase boundaries approach the 
tricritical point. It therefore has the potential of provid- 
ing relevant and well-defined information about the phase 
structure of the NJL model, albeit in a limited region of 
the phase diagram. 

This paper is organized as follows. In Sect. II, we out- 
line the derivation of the GL effective action. In Sect. Ill, 



As in previous works we use the language of the rela- 
tivistic HF approximation to implement the large N limit 
[1^ . The single particle Dirac Hamiltonian with scalar 
and pseudoscalar mean fields S, P reads 



H ^ -7^i9:r + S{x)f + P(x)i7\ 
where the 7-matrices can be chosen as 



(5) 



7° = cri, 7^ = -io'2, 7^ = CT3. (6) 

In HF approximation the grand canonical potential con- 
sists of a single particle contribution plus a double count- 
ing correction, 

{S - mo)2 + P2 



c\x- 



2Ng^ 



(7) 



and has to be minimized with respect to S and P. The 
complicated part is the dependence of the spectral den- 
sity (j{E) on S and P. The spectral density 

(t{E) = TrS{H - E) = -lmR{E + ie) (8) 

TT 

is closely related to the resolvent of H, 



1 H + z 
R[z) = Tr— = Tr— ^ -. 

H — z H'^ — z^ 



(9) 



If the mean fields are weak and slowly varying, we can 
apply the derivative expansion [lH to derive an approx- 
imate, closed expression for a{E) and hence for . It 
will involve a polynomial in S, P and their derivatives 
and may be viewed as a microscopic GL theory for the 
NJL model near the tricritical point. The calculation is 
straightforward and has been done along the lines ex- 
plained in more detail in [l^ [13] . We decompose H and 
as follows, 

H = K + I, ^Hl + V. (10) 

The resolvent can then be expanded formally as 

i?(z) = Tr(X + / + z) (Go£(-1^Go)") (11) 



3 



with 



Go — 



(12) 



We identify the unperturbed part K oi H with the free 
massless Dirac Hamiltonian and similarly Hq with 
and combine the mean fields 5, P into one complex field 

$ = 5 - iP. (13) 

This specifies all operators entering the derivative expan- 
sion (fTT|) as follows, 



K 



-\d^ 

19:^ 



Hi - -dl, 



V = 



m 



$ 

$* 

4$' 



i($')* 1*1^ 



(14) 



As in a previous application of this formalism to baryons 
[l3l all infinities can be controlled with the help of the 
vacuum gap equation 



(1-mo) 1 
— l^TT- = -1^^ 



(15) 



and the following relation for the chiral symmetry break- 
ing parameter 



7 = mo In A, 



(16) 



where we use units where the vacuum fermion mass is 
m = 1 for all 7. Counting the field 4> and each derivative 
dx formally as being of order e, we evaluate all terms 
in the expansion (fTTj) up to order e^. With the help of 
the momentum space technique developed in ,17|] , we find 
after a lengthy calculation the effective action (equivalent 
to the grand canonical potential density) 

^f^g = ao + ai - 2Re$) -f q;2|^'P -t- Q;3lm$($')* 
+ a4(|$|* + |$f )-)-a5lm($" - 3|$|2$)($')* (17) 
+ ae (2|$|6 + 8|$H$f + 2Re ($')'(^*)' + I*'?) 

The coefficients are (7, /i, T)-dependent functions 

an = 1 

" 6 27r 47r 

ai 



1_ 
27r 
1 



"2 = — [ln(47rT) + ReV(2)] 
zvr 



1 



as = 

ai = 

"5 = 

ae = 



237r2T 
1 



26^37^2 
1 

2837r4r3 
1 



Im?/; (1, z) 
ReV'(2, z) 



Im-0 (3, z) 



2i237r5T4 



ReV'(4,z) 



(18) 



expressed through standard di- and polygamma functions 
of 



1 



27rT' 



(19) 



For real Eqs. (jlTHlOp reproduce the effective action 
for the GN model which was derived in [l3[ from the full 
solution. The novel terms (~ as, as) can be tested in the 
chiral limit where the exact <i> is known from the chiral 
spiral 14]. The derivative expansion is particularly useful 
in a case like the NJL model where the exact answer is 
not known. 

In the GN model, it was imperative to go to 6th or- 
der since the terms of order and e'^ both vanish at the 
tricritical point. In the NJL model, this is not the case 
since the tricritical point is now located elsewhere. Here 
it is sufficient to keep the terms up to 4th order (ao...a4) 
for a leading order calculation. The original HF problem 
is then reduced to a complex theory which can be 
treated classically in the large N limit (<& is the HF poten- 
tial, a classical quantity). The Euler-Lagrange equation 
then becomes the (inhomogeneous) complex non-linear 
Schrodinger equation 

a4$" - ia3$' - (ai + a2 + 2a4|$p)$ + ai = 0. (20) 

The solution of this differential equation yields directly 
the HF potential, without need to solve the Dirac-HF 
equation self-consistently. This drastic simplification 
only occurs in regions of the phase diagram where the 
mean field is both weak and smooth. Unfortunately, it 
does not seem possible to integrate Eq. (PU)) in closed an- 
alytical form. Eventually, we shall solve it numerically in 
order to minimize the effective action. Before turning to 
the results, we describe a semi-quantitative variational 
calculation. It has all the salient features of the full nu- 
merical solution while being analytically tractable. 



III. A SIMPLE VARIATIONAL CALCULATION 

In the present section we minimize the (4th order) GL 
effective action 



^'(,ff = ao -I- ai - 2Re$) -I- a2|$P -I- a3lm<I>(<I>') 

|4 , |^/|2\ 



2 _ OT}„ I l,T>|2 

-fa4 (|<i>|^ + |$T) (21) 

using the following variational ansatz 

$ = M-fyle"5^, (22) 

with three real parameters M, A, Q. The ansatz is mo- 
tivated by the fact that it contains the homogeneous so- 
lution {A — 0) and the chiral spiral {M = 0) as special 
cases and by the desire to perform the whole calcula- 
tion analytically. Inserting Eq. ((22|) into Eq. pT|) and 
averaging over one period L = 2tt/Q yields the effective 
potential 

= + ai{M^ + A^ - 2M) + a2{M^ + A^) 
- a^A^Q + a4 [M* + A* + A^{AM^ + Q'')] (23) 



4 



This 4th order polynomial has to be minimized with re- 
spect to M, A, Q. We first vary with respect to Q, 



ag^-cff = ^2 {2a4Q ~ as) = 0. 



(24) 



This condition yields either A ~ 0, i.e., the homogeneous 
phase 

^'hom = ao + ai{M^ - 2M) + 03 Af^ + a^M^, (25) 
or it determines the wavenumber Q to be 

(26) 



In the latter case, we arrive at the effective potential 
governing the crystal phase, 



cryst 



ao 



A' 



2M) + a2{M^ + A^) 

2 

A'^ + 'iAHd^)- -^A^, (27) 
4a4 



which we minimize with respect to A^ next 



4a4 



or 



2a4 



A real value of A is only obtained if 



16al 



ai + a2 
4a4 



Sal 



(28) 



(29) 



(30) 



Keeping in mind this upper bound for M^, we eliminate 
A^ from Eq. (|27p and get an effective potential of the 
single variable M for the crystal phase. 



cryst 



(ai + a2)^ a|(ai + a2) 



a? 



4a4 

2aiAf — ( ai + a2 



8al 



64at 

^]m^ - SaiM^ (31) 
2a4 / 



as 



It is now possible to combine the information about both 
phases contained in Eqs. (|25l3ip into a single effective 
potential of M as follows. For AP > M^^^, '^cS is unique 
and given by *hom, Eq. 1^. For < M^^^, is 
double valued, equalling either ^'hom or V&cryst- Actually, 
in this region 



cryst 



horn 



-4a4(M2 - Ml,^f < 0. 



(32) 



Hence the homogeneous branch can simply be disre- 
garded in the region AP < M^^^ when minimizing 4" off. 
The following effective potential thus contains all the per- 
tinent information about the phase structure. 



-Af2)*„y3t. (33) 



homogeneous 




FIG. 2: Phase boundary at 7 = 0.0005 according to the vari- 
ational calculation based on the ansatz (I22|) . separating the 
homogeneous phase from the crystal phase. A 1st order line 
to the left of the tricritical point Pt goes over into a 2nd order 
line to the right of Pt. The points labelled (a,b,c) refer to the 
effective potentials shown in Fig. 3. 



Although the definition is piecewise, the crystal and ho- 
mogeneous effective potentials and their first derivatives 
with respect to Af match a.t M = zbAfmax, so that our 
final effective potential is a smooth function of AI (the 
2nd derivative is discontinuous at ztAfmax)- This enables 
us to analyse the phase structure in a simple way. The 
system is in the crystal or the homogeneous phase de- 
pending of whether the global minimum of "^cff is in the 
region A4'^ > Af^^x or ^'^^ < -^^max- ^ first order phase 
transition can be identified by an effective potential hav- 
ing two minima of equal depth, one on each side of Afmax- 
If a 2nd order phase transition happens, there is a unique 
minimum which crosses the boundary AI = A^max- Since 
^'cff is a 4th order polynomial in AI, the minima can be 
found analytically by solving a cubic equation. 

We now turn to the results. In Fig. 2, we show the 
phase boundary obtained at 7 = 0.0005, restricting the 
(/i, T) region to a neighborhood of the tricritical point. 
This plot is representative for all 7 <C 1. The most strik- 
ing result is the appearance of a tricritical point Pt sep- 
arating a first order transition line (to the left) from a 
2nd order line (to the right). To illustrate the way in 
which the order was identified, we show in Fig. 3 the ef- 
fective potentials belonging to the points in Fig. 2 labeled 
by (a,b,c) and crossing the phase boundary at three dif- 
ferent temperatures. As explained above, the effective 
potential contains crystal and homogeneous solution to 
the left and right of A/max- One clearly sees the differ- 
ence between the first order transition (a) and the 2nd 
order transition (c). At the tricritical point (b), the 2nd 
derivative ^"g (Af) to the left (i.e., in the crystal phase) 
vanishes in addition to 1st derivative (remember that ^'"g 
is discontinuous). Inspection of the effective potentials 
reveals that the instability can be described as follows: 
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FIG. 3: Evolution of effective potential when crossing the 
phase boundary in Fig. 2 at three different temperatures, illus- 
trating the difference between 1st and 2nd order transitions. 
A minimum to the left of Mmax refers to the crystal phase, 
a minimum to the right of Mmax to the homogeneous phase. 
Figures 3a, 3b, 3c correspond to the points a,b,c marked in 
Fig. 2 

Above Tt when moving from left to right through the 
phase boundary, there is a perturbative instabihty to- 
wards modulating the mean field, 

$ = Af ^ $ = Af -f Ae'^'^ (34) 

M is continuous and A vanishes at the phase boundary. 
This is the standard perturbative scenario also observed 
in the GN model in the upper sheet of the phase diagram 
p^ . Below the mean field changes discontinuously 

$ = M $ = M' + Ae"5^ (35) 



FIG. 4; Two-dimensional phase boundary in (7, ^, T) space, 
separating a massive Fermi gas from the chiral crystal, ac- 
cording to the variational calculation. Notice that only the 
vicinity of the tricritical point is shown. 



where M ^ M' and A jumps from to a finite value. 
This is remarkably different from the GN model which 
shows another type of instability here, namely towards 
formation of a single baryon in a 2nd order transition 
[l^ . Repeating this calculation at several values of 7 
finally yields the phase diagram shown in Fig. 4. It has 
two sheets of 1st and 2nd order transitions separated by 
a tricritical line and nicely matches the phase diagram of 
Fig. 1 in the limit 7 — s- 0. 

Let us now try to understand how exactly the tricritical 
line approaches the tricritical point = (0,0, Tc)- On 
general grounds one expects some power law behavior in 
variables which vanish at the phase transition. We there- 
fore first replace the temperature by the shifted variable 



T = ^T, - T. (36) 



In Figs. 5a and 5b, we show the phase boundaries for a 
sequence of 7-values in the (^, T)- plane (a projection of 
the surface in Fig. 4 onto the (/x, T)-plane) and in the 
(/i, r)-plane. This is not yet very illuminating. However, 
if we now rescale both axes in Fig. 5b by a factor 7"^/"^, 
we obtain the striking result shown in Fig. 5c. We see 
that this particular rescaling removes the entire 7 depen- 
dence, at least within the resolution of our plot. Thus 
the information contained in the 3d plot Fig. 4 could 
have been inferred from a computation at a single value 
of 7, e.g. the curve shown in Fig. 2. The tricritical 
line is just one point in the scaling diagram located at 
(Ai7"^/^r7-l/3) « (1.488,0.488). In the next section, 
we show that this scaling behavior has nothing to do with 
the specific variational ansatz, but can be derived more 
generally from the effective action and simplifies the full 
calculation considerably. 
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the coefficients a;, Eq. (jlSp . around the tricritical point, 



0.6 0.8 



-1/3 



FIG. 5: a) Projection of the 3d phase boundary of Fig. 4 onto 
the (/i, T) plane, b) Same plot, replacing the variable T by 
T — VTc — T. c) Rescaling both axes by a common factor 



-1/3 



removes the 7-dependence. 



IV. TRICRITICAL BEHAVIOR IN THE FULL 
CALCULATION 



We first explain the reason behind the simple scaling 
law observed in the above variational calculation. We 
restrict ourselves to values of the independent variables 
(7,/z,t) of the order 



7 



(37) 



This will be justified afterwards by showing that the tri- 
critical line indeed lies in this region. We then expand 



7 
7 

«3 = -^C(3)e-2^/i + 0(e3) 



,-2C ,,2 



-C2 



+ O(e^) 



8tt 



-2C 



(38) 



Returning to the Euler-Lagrange equation, Eq. (pn|) . di- 
viding by Q!4 and using the expansion ([35]), we find to 
leading order 



- 4i/i$' - 4/i 



Tc 



7-2 + 2|$n $ + a7 = (39) 



where we have introduced the transcendental constant 
16e2C 



6.032 



(40) 



7C(3) 

We next define a rescaled scalar field 

$(x) =7^/3^(71/3^) (41) 

with 

$'(x) = 72/3<^(7i/3x), $"(x) = -fifii-f^/^x), (42) 

where ip denotes the derivative of (p with respect to the ar- 
gument ^ = 7I/3X. Inserting this ansatz into the differen- 
tial equation (|39p and employing appropriately rescaled 
variables for chemical potential and temperature (the 
constant factors of 0(1) are introduced to simplify the 
notation below). 



v = 27" 



-1/3 



a=J-7-^/^7 



the 7-dependence drops out, 

(p — 2\VLp — (l/^ — cr^ 



2\ipY) ip + a = 0. 



(43) 



(44) 



This confirms that the 7-dependence in the original equa- 
tion was spurious and can be trivially reconstructed from 
a calculation at a single value of 7. Under these transfor- 
mations, the grand canonical potential scales as follows. 



off 



+ (45) 



7 
27ra 

-\-2v Ira ip{ip) 



The value of 7 only enters as an overall factor. As a con- 
sequence, when drawing the phase boundary in the (z/, a) 
diagram rather than the (/x, T) diagram, one expects a 
universal curve. Needless to say, this only gives the lead- 
ing asymptotic behavior of how the tricritical point is ap- 
proached, since we truncated the Taylor expansion of the 



coefficients of the GL effective action around this point 
after the leading terms. 

We now turn to the fuU computation of the universal 
phase boundary which replaces Fig. 5c from the varia- 
tional approach. The calculation is done differently for 
the 2nd and 1st order phase transition. The 2nd order 
phase transition can be treated perturbatively. In this 
case, we set 

$(a;) = M + ip{x) (46) 

and keep only terms up to 2nd order in 1^9(2;) when com- 
puting the grand canonical potential. Denoting the dis- 
crete Fourier components of (p{x) by ipi, the grand canon- 
ical potential averaged over one period will assume the 
form 

/ dx'i'^B = C + ^Fu,^i^i,. (47) 

Vanishing of det F signals the perturbative instability. F 
being block diagonal (only the are non-zero due 

to spatial averaging), it is sufficient to consider a single 
Fourier mode in order to detect the instability. Since the 
period is also a variational parameter, we may further- 
more assume \i\ = 1 without loss of generality. As a 
result, the perturbative calculation is similar to the toy 
model of the previous section, except that one has to use 
the more general ansatz where the field $ traces out an 
ellipse rather than a circle in the complex plane, 

$(x) = M + Acos{Qx) + iBsiniQx). (48) 

We expand the effective action to 2nd order in A, B and 
study where the homogeneous system gets unstable. The 
information which we get in this manner is the boundary 
of this perturbative instability, the ratio of the main axes 
of the ellipse, R — A/B, and the wavenumber Q in the 
unstable direction at the boundary. Let us introduce 
once more rescaled parameters for the mass M and the 
wavenumber Q, 

M = m-f^/^, Q = 97^/3. (49) 

We then find the following parametric representation of 
the critical curve, 

2miy^ = a + 2m^ + ^/a(a + Am^), 

2ma^ = -a + 6m^ + ^y a(a + Am^). (50) 

The parameters q, R of the unstable mode are given by 

2 2,2^2 
q = u + a — 4to , 

„ + 0-2 — Am? 

R = ^ 5 . (51) 

+ 

The curve (I50p is shown in Fig. 6. This perturbative 
calculation does not teach us where the tricritical point 
is, and therefore where the 2nd order phase boundary 
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1 4 6 8 10 12 14 16 18 



FIG. 6: Perturbative phase boundary in the rescaled variables 
(v, o) defined in Eq. (|43[) . computed from the parametric rep- 
resentation ifSU)) . 



1.2 



0.8 

R 

0.6 




0.4 
0.2 



" 2 4 6 8 10 12 14 16 18 

V 

FIG. 7: Ratio R = A/B oi the main axes of the ellipse traced 
out by $ in the complex plane, right at the perturbative phase 
boundary. 



ends (it starts at small cr, large i.e., from the bot- 
tom right of the plot). In Fig. 7 we show the ratio R 
characterizing the eccentricity of the ellipse. In the re- 
gion of interest later on, it deviates from 1 (the value 
assumed in the variational calculation) by less than 6%. 
The reduced wavenumber q is equal to v for all practical 
purposes so that we refrain from showing a plot. Going 
back to the original variables, this perfectly agrees with 
the wavenumber Q — 2fi oi the chiral spiral at 7 = 0. 

We now turn to the first order transition line and the 
tricritical point. Here, perturbation theory is of no help 
since the transition is discontinuous. We therefore have 
to integrate the rescaled Euler-Lagrange equation numer- 
ically. In order to explain our method, it is useful to rein- 
terpret the problem as one of classical mechanics. If we 
identify the spatial coordinate ^ with time and S, P, the 
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real and imaginary parts of $, with cartesian coordinates 
in a plane, the effective action becomes the Lagrange 
function of a (planar) non-relativistic particle moving in 
a static electric and magnetic field, with the equation of 
motion 

m?=q{E + fxB^. (52) 

The electrostatic potential V{r) {E = —VV{f)) is given 
by 

^V{f) = ax-^{u'^ -a^y -^r"^. (53) 

Depending on the parameters {i^, a) it looks like a tilted 
inverted Mexican hat (<t > z^) or a tilted sugarloaf 
((T < v). The inversion is due to the interchange of 
space and time, chiral symmetry is rotational symme- 
try in the plane, and the term ax arises from explicit 
chiral symmetry breaking. The magnetic field is homo- 
geneous and perpendicular to the plane of motion with 
magnitude proportional to the chemical potential, 

—B = -2iye,. (54) 
m 

There are of course infinitely many solutions of Newton's 
equation, depending on the initial conditions. The ho- 
mogeneous solution for (p {f = m) corresponds to the 
time independent solution where the particle is sitting 
at the maximum of the potential. Since we expect ip to 
be spatially periodic in the crystal phase, we search for 
periodic orbits of the classical particle and minimize the 
action with respect to the initial conditions. We proceed 
as follows. First we notice that energy 

E = y f 2 ^ ^y(^) (55) 

is conserved. In the original problem, this implies that a 
certain function of (p and (fi is ^-independent, 

(|<^|2 - (z/2 - a2)|<^|2 - |(^|4 + 2aRe^) = 0. (56) 

For symmetry reasons, we can choose without loss of gen- 
erality the initial conditions 

where vo can be computed from the energy. 



vo = J-{E-qV{xo,0)). (58) 
V m 

We then vary Xq and solve the differential equations nu- 
merically until we find a periodic orbit. This requires 
an iterative procedure. The energy E is treated as a 
variational parameter. The calculation is repeated for 
different energies until the minimum of the action has 
been found. 



Im (p 











\ \ \ \ 












-2- 





FIG. 8: Sample periodic orbits of a classical charged parti- 
cle in an electromagnetic field, equivalent to solutions of the 
Euler-Lagrange equations for the rescaled GL action. The 
dot at 95 « 0.78 corresponds to the time independent solu- 
tion right at the tricritical point (i^, a) — (2,99,1.54). The 
curves arc obtained for the parameters (j^, (j) — (2.93,1.62), 
(2.58,2.15), (2.36,2.57), (2.20,2.93), (2.09,3.25), (1.99,3.54)^ 
from inside out. Since the shape of the trajectories is reflec- 
tion symmetric only the lower half is shown. 



We now turn to the results obtained with Maple. We 

find that the energy E needed to minimize the action is 
always very close to the maximum value of the poten- 
tial. Periodic orbits are easily found at small a but are 
increasingly difficult to find at larger a since one needs 
to fine-tune the initial coordinate xq to a very high pre- 
cision. This is due to the fact that the inverted Mex- 
ican hat shape becomes more pronounced and the par- 
ticle travelling along the rim will in general fall inside 
or outside before making a full turn. Some representa- 
tive orbits are shown in Fig. 8. The shape agrees very 
well with the circle assumed in the variational calcula- 
tion. The most conspicuous discrepancy is the behavior 
near the right end for large radii where the circle gets 
deformed. This happens close to the maximum of the 
inverted Mexican hat potential. A closer look at the 
time dependence shows that the particle spends more 
time in this region, in contrast to the assumed uniform 
motion underlying the variational calculation. The tri- 
critical point can be rather easily identified by the disap- 
pearance of the non-perturbative solution. Fig. 9 shows 
the rescaled critical line. We have combined the numer- 
ical result for the 1st order line with the analytical cal- 
culation of the 2nd order line. The rescaled tricritical 
point lies at Ft ~ (2.99, 1.54), very close to the value 
(2.98, 1.59) of the variational calculation. As explained 
above, the curve shown in Fig. 9 is all what is needed 
to reconstruct the 3-dimensional phase diagram in the 
vicinity of the tricritical point, see Fig. 10. This gives a 
first glimpse of the full phase diagram of the massive chi- 
ral GN model which is still largely unknown. Although 
we cannot rule out further phase transitions, we believe 
that it contains already the crucial information about 
the existing phases and the order of the phase transi- 
tions. By our derivation it should correctly describe the 
asymptotics of how the tricritical point is approached as 
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1st order 



2nd order 



FIG. 9; "Universal" phase boundary in rescaled variables 
obtained from perturbative, analytical (2nd order) and non- 
perturbative, numerical (1st order) calculations. 




0.001 "iV 0.02 



FIG. 10: Like Fig. 4, but result from numerical minimization 
of the GL effective action. The whole surface has been recon- 
structed from the single curve in Fig. 9 by using scaling. The 
curves 7 = const, would all have the same shape in a {fi, r) 
plot, but a different scale. The second set of curves drawn 
onto the surface connects points related by scaling. The nar- 
row gap between the phase boundary and the T-axis is an 
artefact of the limited range of a when determining the 1st 
order curve in Fig. 9. 



we let /i ^ 0, 7 ^ 0, T ^ Tc. To illustrate this claim wc 
mention that the tricritical line along which the order of 
the phase transition changes behaves asymptotically (for 
7-^0) like 



1 

2^t7 



1/3 



T = Tr 1 



a 



(59) 



where now all the dynamical information is encoded in 
the pair of numbers k, 1.54, Ct ~ 2.99 which had to be 
determined numerically. 



V. COMPARING MODELS WITH DISCRETE 
AND CONTINUOUS CHIRAL SYMMETRY AND 
CONCLUSIONS 



To put the results of this work into perspective, we 
remind the reader that the GN model with discrete chi- 
ral symmetry and the NJL model with continuous chiral 
symmetry, Eq. ([T]) , were believed to have identical phase 
diagrams for many years, see Fig. 2 of Ref. [ll|. This 
was the result of the unjustified assumption that only 
homogeneous phases would matter. The correct phase 
diagrams were first obtained for the massless models, see 
Fig. 1 of the present work for the NJL model and Fig. 4 
of Ref. [ll| for the GN model. Subsequently the phase 
diagram of the massive GN model has been constructed 
(Fig. 5 of Ref. [ll|). Here for fixed 7, the solitonic crystal 
phase is bounded by two 2nd order critical lines. They 
may be characterized by the following instabilities: If one 
lowers the temperature at sufficiently high chemical po- 
tential, the homogeneous condensate (fermion mass) is 
unstable towards a modulation with a finite wavenumber 
Q but vanishing amplitude. If one raises the chemical po- 
tential at sufficiently low temperature, the system is un- 
stable with respect to formation of a single kink-antikink 
baryon, i.e., a finite amplitude but vanishing wavenum- 
ber. These two lines are joined in a cusp at a tricritical 
point. If one approaches the tricritical point along the 
first phase boundary, Q goes to zero; along the second 
boundary, the amplitude goes to zero. All transitions 
are continuous. 

In the case of the massive NJL model, we do not yet 
have such a complete picture. As a first step we have 
presented here a study in the vicinity of the tricritical 
point, perhaps the most interesting part of the full phase 
diagram. Again we find at fixed 7 two phase boundaries 
meeting at a tricritical point. To the right of this point, 
we have identified a 2nd order critical line with a per- 
turbative instability similar to the one seen in the GN 
model. Here, in addition to a spatial modulation, the 
mean field acquires an oscillating imaginary part from 
the pseudoscalar condensate. Another important differ- 
ence is the fact that the wavenumber Q of the unstable 
mode does not vanish as we approach the tricritical point. 
At the tricritical point, there is no cusp, but the phase 
boundary continues smoothly. In contrast to the GN 
model it now turns into a discontinuous 1st order phase 
transition. The system becomes unstable with respect to 
formation of a chiral crystal with finite period and ampli- 
tude, rather than creation of a single baryon. This is sur- 
prising, since the massive NJL model is known to possess 
solitonic baryons. At T = one would definitely expect 
a 2nd order transition at a chemical potential equal to 
the baryon mass, but apparently the transition becomes 
1st order at finite temperature. 

In conclusion, the study of the NJL thermodynamics 
is not just a mere repetition of that of the GN model, but 
involves interesting new challenges and physics facets. 
Due to the complex order parameter and the fact that 
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one can apparently do much less analytically, it is quite 
non-trivial. Nevertheless it would be interesting to de- 
termine the complete phase diagram. Exactly soluble 
relativistic field theories with a rich phase structure are 
so rare that one should not simply ignore this problem. 
Moreover, there may be applications in condensed mat- 
ter physics like in the case of the GN model. Thus for 
instance, relativistic four-fermion models with a contin- 



uous chiral symmetry have been discussed in the context 
of incommensurate charge-density waves some time ago 
[H, [l^, [13] . Spiral order parameters have an even longer 
history, starting with the classic paper o n sp in density 
waves in the electron gas by Overhauser (2l!|. Last not 
least, let us mention that various chiral crystals have been 
studied in 3-1-1 dimensional quantum chromodynamics, 
see e.g. the recent work [22i] and references therein. 
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